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Abstract: This paper presents a mathematical model for a sequential imperfect preventive maintenance 
policy which considers stochastic factors affecting the failure rate and service life of a repairable 
deteriorating system. The failure rate and maintenance expenses depend on a number of influencing factors, 
such as the age, condition of the equipment, plant capacity utilization, mode of operation as well as the 
applied maintenance concept. With the failure rate following the Weibull distribution, the model investigates 
the combination of the age reduction and failure rate increase coefficients to determine the optimal intervals 
between preventive maintenance schedules and overhauls in order to minimize the expected costs over an 
infinite time horizon. The optimal policies are presented and numerical examples are also given. 
Keywords: Preventive Maintenance (PM), Age Reduction, Hazard rate, Correction Factors, Cost 
Minimization 



I. INTRODUCTION 

The most commonly used maintenance models for the failure process of a deteriorating system are perfect 
maintenance actions which restore the system to 4 as good as new’ and imperfect maintenance which restores a 
system to a state between ‘good as new’ and ‘bad as old [1] .The link between the two concepts is done by corrective 
maintenance (CM) or minimal repair which restores the system to the state it was prior to failure. In practice, the 
reality is between the two extremes of perfect and imperfect maintenance. The perfect maintenance may not yield 
a proper functioning system which is ‘as good as new’ and the imperfect maintenance seems to be too pessimistic 
in repair strategies as it depends on the availability of resources. Timely preventive maintenance actions can inhibit 
further equipment degradation and reduce the frequency of failures thereby extending its service life. Ideal 
maintenance policies aim to optimize system reliability, availability and safety performance at optimum costs 
[2], [3 ],[4],[5]. 

The optimization of maintenance for repairable systems was initiated by [6]. Since then, a number of 
maintenance models and strategies for maintainable systems have appeared in literature. An imperfect preventive 
maintenance model was introduced by [7] which described the age reduction effect and that preventive maintenance 
reduces the failure rate of a system to a fraction of its value just before PM. Then another model by [8] where a 
system undergoes PM with a different failure distribution between preventive maintenance actions was considered. 
An age-dependent preventive maintenance model [9] was presented and it was observed that, the selection of an 
appropriate maintenance strategy directly depend on the previous operating parameters. 

In [10] the problem of optimal PM policies was investigated and assumed that after each repair the 
system’s age becomes zero(as good as new) while the failure rate increases by the increasing of the number of 
repairs carried out. 

A periodic PM model [1 1] for which the maintenance slows the degradation process of a system, while the 
hazard rate is monotone is discussed. Authors in [12] formulated maintenance cost policy where PM is performed 
based on predetermined maximum failure rate, and failures are corrected by minimal repair. 

Improvement factors in the hazard-rate function were adopted in [13] and analysed two corresponding 
imperfect PM models: PM reduces the hazard rate while it increases with the number of PM actions and that PM 
reduces the age. In [1] are listed methods and techniques used for modelling imperfect preventive maintenance and 
observed that, maintenance models would be valuable if they are capable of incorporating all factors that affect a 
repairable system. 
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A hybrid PM model comprising a combination of the hazard rate and age was proposed in [14]. They 
further, [15], studied the problem of optimizing PM for repairable systems with two categories of failure. A similar 
work was proposed by [16], they developed an age-based hybrid model for imperfect PM with the two categories of 
failure modes, i.e. maintainable and non-maintainable failure modes. An optimal periodic PM policy for the 
deteriorating system based on [1 1] was proposed in [17], i.e, a maintenance model where PM slows down system 
degradation while the failure rate remains monotonically increasing was analysed. Optimal maintenance options 
[18] for continuously monitored multi-component systems based on Markov deteriorating processes were studied. 
In [19] the author considered the imperfect effect of PM activities and gave an availability-centred preventive 
maintenance model for multi-unit systems, which is based on sequential PM theory. 

Several preventive maintenance models that considered the effect of imperfect maintenance on effective 
age of a repairable system are presented in [20]. And [21] further built a model where each PM action reduces the 
effective age but not restoring it to an 4 as good as new’ state. A sequential imperfect preventive maintenance policy 
with random maintenance quality under reliability limit was proposed by [22]. And [23] proposed improvement 
factors in a sequential maintenance model. Similar work was done by [24] while considering phasic sequential PM. 
A stochastic maintenance policy for a degrading system over a finite period was considered by [25]. 

The general form of maintenance optimization analysed is either based on cost, age reduction, and failure 
rate or reliability/availability analysis. This paper looks at a specific function using both the concept of increasing 
failure rate and age reduction for a deteriorating maintainable system to determine intervals between maintenance 
actions in order to minimize the related maintenance costs. 

II. MODEL DESCRIPTION 

In [14] adopted correction factors for failure rate, h(t), and time interval between preventive maintenance 
actions, i.e., the failure rate during the subsequent period is ah(t), where h(t) was the failure rate prior to the last 
maintenance action, a > 1 is a correction factor and t>0 is the time after the previous maintenance action which is 
reset to zero at each PM. Time interval between two maintenance actions is reduced to bt after each maintenance 
action, where b < 1 is the reduction factor. According to [13], for a failure rate h(t), t E(0, fi), the preventive 
maintenance action at t\ determines a new failure rate g(t ), t E(7 1? t 2 ), which depends on the previous failure rate and 
the quality of maintenance action, i.e. failure rate reduction. The proposed form of g(t) uses both the concept of 
increasing failure rate and reduction of preventive maintenance interval according to system usage. The failure rate 
model assumes that the failure rate immediately after preventive maintenance action increases more quickly than it 
did before the PM action. The age model assumes that there is an effective age reduction immediately after the PM 
action, and the failure rate continues to be a function of the effective age. Based on the above observations, the 
model captures both by how much the effective age is reduced just after PM action and how much faster the failure 
rate function increases when the system is back into operation. The model failure rate function after the first PM is 
written as: 

-f x) = ah(bt L + x) (1) 

Where x E(0,t2-ti). For a - 1 the model is resized to the one of life reduction, and b = 0 to increasing the failure 
rate. 

Using the failure rate in (1), optimal preventive maintenance policies are developed, with a major implication, i.e. 
giving up the classic, non-economic approach of constant intervals of preventive maintenance execution in order to 
avoid costly sudden failures, with the following assumptions. 

a. The repairable system under consideration deteriorates with increased usage and age and is regarded as a 
single-unit system without consideration of dependencies between different components. The planning 
horizon is infinite. 

b. The equipment undergoes minimal repair upon failures between two adjacent PM intervals. The times for PM 
and minimal repair are negligible. 

c. The imperfect PM is performed at a sequence of time fi, t 2 , .... t N _i and is overhauled at t N . 

d. The failure rate of the system is increasing despite the PM actions. 

In the proposed maintenance model, the following notations will be used: 

-hit) is the failure rate and H(t) is cumulative failure rate of the system, 

-X is the maximum value for failure rate; 

-x- K , k = 1,2, iV, length of scheduled intervals between PM actions, 

-yk is the effective age of the system just before r: r * preventive maintenance action, 
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-N is the number of preventive maintenance actions, (. N th action is perfect maintenance which is overhaul) 

-Q-k is the failure rate adjustment factor after k th preventive maintenance action with 
1 = < G-i < a 2 < ... < o N _ L ; cine A- k = IIi=o , k = 1,2, .,..,A r 

-b k is effective age adjustment factor after rc f ' : preventive maintenance action, 

0 = b 0 < h 1 < b 2 < — < b N _ ± < 1; 

-c m is the corrective maintenance cost; c tl is preventive maintenance cost; c r is overhaul cost; C is average 
cost rate. 

It is assumed that the system is subjected to preventive maintenance at times t u t 2 ... t N -i and overhauled or 
replaced at t N _ Corrective maintenance is executed if failures appear between adjacent preventive maintenance 
actions. The overhaul at the time t N doesn’t restore the system to ‘as good as new’ state. The system has a failure 
rate A k h(t) between the number t k .j and t k preventive maintenance actions, i.e. in the range (t k _ i, t k ). The effective age 
of the system is b k _iy k _i immediately after the k- 1 preventive maintenance action, and it becomes y k = x k + b k _ pc k _i 
+ ...+ b k _ib k _ 2 ...b 2 biXi after the number k preventive maintenance i.e. the age changes from b k _cy k _ i to y k in the range 
(t k . i, t k ). 



y k = x k -v b k .\y k .\ 
x k = y k - b k _iy k _i 

Given that 



xi =yi 
x 2 = y2~ bfy ! 
x 3 = y3~ b 2 y 2 
x n -yN~ bui yN-i 

Then 

Cycle = Si = 21=1 Cl - iftVk -I- y N 



Average cost 

C(y v y 2 ...y K } 



C r +Ct;( j V + c mEij=i (j'ft.) ^C^fe-lJ'fe-l.' 1 ] 



( 2 ) 



(3) 



2.1. Model A - Cost Minimization with Failure Rate limit 

The objective of this model is to perform preventive maintenance action before the failure rate, X, exceeds 
a prescribed value, i.e. at the time t k (k = 1, 2, ..,N). As the system failure rate is strictly increasing, preventive 
maintenance intervals will also keep changing to avoid system sudden failure. 

That is, 

A = A* k = 1,2, . . , N (4) 



Solving equation (4), we get y k (k = 1, 2, ..,N) as a function of A. Substituting (4) into (3) and differentiating with 
respect to A = 0, (t. e. = 0 ] deads to 



Cvt J - * b k 7il> . Kv n ) 

yjV-l r 1 = 

k=i-A k h(j k ) Axhirx) 



C 

c m 



(5) 



Solving (5) with respect to X and get X as a function of N, we find N which minimizes the left hand side of equation 

(5). 

Based on the above, the preventive maintenance interval algorithm is derived as follows: 

1. Solve equation (4) and find y k with respect to 2; 

2. Substitute solutions in step 1 into equation (5) and determine A ; 

3. Choose A to minimize the left hand side of (5), where y k ( k = 1, 2,..,N) is from step 1 & 2; 

4. Compute y k (k= 1,2 ... N) from steps 1 and 2, for the value of N from step 3; 

5. Compute x k = y k - b k _ x y k _ u k = 1, 2... N. 
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2.2. Model B - Cost Rate Minimization with Age Reduction 

In the age-reduction model, it is assumed that the effective age of the system after the k th PM is reduced to 
b k y k if it was y k before the k th PM as the number of PM actions on the system increases. The decision variables are 
N and the maintenance ages y h y 2 ... yN that minimizes C(y h y 2 ... y n) in equation (3). 

3 C 

differentiating equation (3), with respect to y k — = 0, we get 

’ 



- ^jc+ifcjcAC&jtyfe) = A w (l - k = 1,2 N - 1 (6) 

And c m A N h(y N ) = C (7) 

Substituting each solution y k of equation (6) into (7), we get 

Efctfl -b k )y k + y w ] - YLl A k (8) 

Where each y k (k = 1, 2... N - 1) is a function of y N . 

The algorithm to find the optimal maintenance intervals is the same as in the failure rate limit model described 
earlier. 



III. NUMERICAL EXAMPLES 

The repairable system under consideration is an Ore Grinding Mill which is subject to preventive 
maintenance actions every 630 operating hours and follows the Weibull life distribution whose scale and shape 
parameters are, p = 6,148xl0‘ 9 and a = 2,462 respectively. 



hit) = 

3.1 For model A, 

solving the failure rate equation (4) with respect to y k , we have the solutions as expressions of X: 

» = * = 1 ' 2 '--' w < 9 > 

Substituting the values of y k in (9) into equation (5) and solve it with respect to X, we get the values of X with respect 
to N. 






C T -HTpCjV- D 

(T - ^ ^9 



Where E(N ) = - a k fafiT 3 " 1 + A/' 1 

The left hand side of equation (5) becomes 



_i_ 

B a 



Ctt +£ p(l¥- 0 



lEiN) la 
fM 



Where 



F Gv) = Sf-JCl - i>v) V" 1 + V -1 



( 10 ) 



( 11 ) 



Minimization of equation (1 1) is equivalent to minimization of: 

s-1 i 
[ff r +Cp(jV-D] * [EMJa 

fM 

The optimal value N* for N is found by the following inequalities: 




Q(N+ 1) > Q(N) and Q(A0 < Q(N~ 1) 



This implies: 
Where 



11/OVD > - and IVOV - l) < - 



wOv) = 



\£[N + D]a- i i 



i a i a 

[EM] [F( j¥+ - \E (jY+ [f ( AO] 



- 0v - 1 ) 
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3.2 For model B 

According to the described model, equation (6) becomes: 

-^ + i MC&itffc)"- 1 = Vi - - fc = i,z. JV - 1 

Solving the above equation we get 

"- 1 

C T +Op(iV- iJ S 



Substituting (12) into (7) leads to 



y ? N = 



| L A % 1 = 



where d k = [ (L ^ 1" 1 



(12) 



(13) 



Then/l w A(y w ) becomes: 

P 



V [ c 



I C r + c v Qv - 1) Z 



{*^( 1 -;)[' I * +A r 1 






J-l 

Hi 



Minimization of ¥ is equivalent with minimization of the function 



B(N') = 



c r 4- CV - 1) 



l 

A*- 1 



N = 1,2 



j-W-i d 

■jV Ti Ji=i u Ji 
Optimal value A* for A is found solving the inequalities: 



The two inequalities involve: 



B(N +1) > B(N) and B(N) <B(N~1). 



DW>— and D CA - l] < — 



(14) 



D n = ^ 



4 tf-i i W-ij 



A S ~ 1 — 4 iI ~ 1 4- d 

jy+i t a N 

The optimal preventive maintenance interval has the following form x k = y k — b- K _ l y k _ = 1,2, ...,A 
where is given in (12) and (13). For accurate results of x k , (k = 1,2... A), the economic parameters of c m , c p and 
c r are taken from the repairable system operation data. For analytical purposes, the ratios c/c p and c/c m are used. 

For this analysis we use the coefficients a k = (6k+l )/(5k-\-l ); b k = k/(2k+l ), c/c p = 2, 5,10,20,50 ; 

c/c m = 4; recommended in [14]. 

The two models show the reduction of maintenance intervals between two successive actions which 
consequently reduces the average cost in equation (3). Table 1 and Table 2 gives the optimal preventive 
maintenance schedules for both models and same input parameters are considered in both cases. 
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Table 1. Optimal PM Schedules for Model A 



Cr/Cp 


2 


5 


10 


20 


50 


N 


1 


2 


3 


4 


5 


Scheduled Preventive Maintenance Intervals 
(hours) 


Xl 


1663.3 


2292.6 


2619.5 


3089.4 


3991.3 


x 2 




1231.7 


1484.2 


1750.4 


2261.5 


X 3 




1006.3 


1159.9 


1367.9 


1767.3 


x 4 




893.5 


968.4 


1142.1 


1475.6 


X 5 






828.5 


977.1 


1262.4 


*6 






717.05 


845.7 


1092.6 


X 7 






620.6 


736.6 


951.7 


Xg 








644.1 


832.1 


x 9 








564.6 


729.4 


XlO 








495.7 


640.4 


Xll 










562.9 















Table 2. Optimal PM Schedules for Model B 



Cr/Cp 


2 


5 


10 


20 


50 


N 


1 


2 


3 


4 


5 


Scheduled Preventive Maintenance Intervals 
(hours) 


Xl 


1520 


2599.2 


2951.1 


3839.9 


5421.1 


X 2 




1386.7 


1575.1 


2049.4 


2893.3 


X 3 




1079.1 


1225.8 


1594.9 


2259.2 


x 4 




1466.1 


1020.4 


1327.7 


1874.4 


X 5 






870.9 


1133.2 


1599.9 


X 6 






752.4 


979.0 


1382.1 


X 7 






1118.4 


684.5 


1202.1 


Xg 








597.8 


1049.8 


X 9 








522.1 


919.3 


Xio 








798.04 


806.5 


Xll 










708.5 


Xl2 










623.02 


Xl 3 










548.3 


Xl4 










482.8 


Xl5 










752.4 



Table 1, shows the reduction of maintenance intervals with increase of the system failure rate in order to 
avoid sudden failure. This shows that, the studied equipment is subject to degradation during its operation, and the 
preventive maintenance, minimal repair and subsequent overhaul actions do not restore the system to 4 as good as 
new’ state. The failure rate and maintenance costs should be analysed at the beginning of each operation cycle 
depending on prevailing operating conditions. 

From Table 2, it is observed that, the scheduled preventive maintenance intervals decreases for all values 
of x k except the latest time element which shows an increase. The system reaches an optimum age beyond which the 
cost of preventive maintenance will not be economical as the intervals become shorter and the maintenance cost 
start increasing. This coincides with [13], “It would be reasonable to do frequent preventive maintenance with age, 
and it is recommended to do the last preventive maintenance as late as possible because the system should be 
overhauled or replaced at the last PM”. This implies that, during the last operation interval, equipment can be 
utilized as long as it satisfies the minimum operating requirements as it awaits overhaul or replacement. 
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IV. CONCLUSIONS 

This paper presents the ways to develop an appropriate preventive maintenance policy for a deteriorating 
repairable system. Two models, which minimize the expected total cost, for an Ore Grinding Mill whose failure 
rate is according to Weibull distribution law, were studied. The analysis combined the effect of the age reduction 
and the failure rate increase coefficients. Therefore, maintenance actions do not only reduce equipment life but also 
change the failure rate distribution. The proposed approaches provide maintenance managers with decision making 
techniques in the light of all available relevant information and data about the system deterioration. The outcome of 
the case study doesn’t show significant variations between the two models. However, use of the failure rate limit 
model can ensure that reliability of the equipment is kept above a predetermined level to avoid the risk of sudden 
functional failure so that the costly repair and downtime expenses are minimized. 
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